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AN ANALYTICAL STUDY OF THE EFFECTS OF JETS LOCATED 
MORE THAN ONE JET DIAMETER ABOVE AWING 
AT SUBSONIC SPEEDS 

By Lawrence E. Putnam 
Langley Research Center 

SUMMARY 

A procedure has been developed for calculating the effects of blowing two jets over 
a swept tapered wing at low subsonic speeds. The algorithm used is based on a vortex- 
lattice representation of the wing lifting surface and a line sink-source distribution to 
simulate the effects of the jet exhaust on the wing lift and drag. The method is limited 
to those cases in which the jet exhaust does not intersect or wash the wing. The pre- 
dictions of this relatively simple procedure are in remarkably good agreement with 
experimentally measured interference lift and interference induced drag. 

The results of the analytical study have indicated that substantial increases in 
wing lift and reductions in wing induced drag can occur when jets are blown over a wing. 
The magnitude of these effects increases with an increasing ratio of jet-exit velocity to 
free -stream velocity. The results also indicate that the interference effects are the 
most favorable when the nozzle exits are located ahead of or at the wing leading edge and 
inboard near the wing center line. The present method does not correctly predict the 
effects of the vertical location of the jets at jet center-line heights of less than approxi- 
mately 1.5 jet-exit diameters; however, as the jet vertical location is reduced to this 
height, the favorable effects of blowing jets over a wing increase. 

INTRODUCTION 

Recently, there has been renewed interest in locating the propulsion systems of 
future jet transport airplanes forward and above the wing. This unconventional engine 
location can take two forms. In one form, the jet is exhausted so that it attaches to the 
upper surface of the wing to take advantage of Coanda turning of the jet exhaust as it 
passes over the wing -flap system to achieve high lift. (See refs. 1 to 4.) In the second 
form, the jet exhaust is raised so that the exhaust flow does not attach to the wing in 
order to avoid the cruise drag penalties associated with the first form. In both cases, 
the jet noise propagated below the airplane during take-off and landing may be reduced 



by the shielding effects of the wing. (See refs. 5 and 6.) To date, however, there is 
very little experimental information of relevance, except for references 7 to 9. Further- 
more, there are only a limited number of theoretical methods (refs. 10 and 11, for 
example) available for use in assessing the effects of jets located above a wing on the 
subsonic cruise and climb performance of future jet transport airplanes. 

The purpose of the present study was, therefore, to develop an analytical method 
by which the effects of jets located above a wing on subsonic cruise and climb perform- 
ance could be investigated. An elementary analytical procedure has been developed to 
calculate the effects of blowing one or two jets over a swept tapered wing. This method 
is based on the vortex -lattice theory for the wing and on the theory of reference 12 for 
a jet exhausting into a subsonic stream. As a result of the basic assumptions, the 
method is limited to the case in which the jet exhaust does not intersect or wash the 
wing. As long as the jet is sufficiently high above the wing so that it does not violate 
this condition, the effects of the following on the lift and drag induced by blowing jets 
over a wing at subsonic speeds can be analytically investigated: ratio of jet to free- 
stream velocity, location of the jet exhaust relative to the wing, nozzle -exit diameter, 
wing leading-edge sweep, wing aspect ratio, and wing taper ratio. 

SYMBOLS 


A e exit area per nozzle 

b wing span 

axial -force coefficient 
Cj) i induced-drag coefficient 

= (^D,i)j e t on " ( CD >i)jet off 
lift coefficient 

- (^L)j e t on “ (C-L)j e ^ 0 ff 

C N normal-force coefficient 

c wing chord or chord of elemental vortex -lattice panel 

c av average wing chord, S re f/b 
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c e wing chord at rj Q 

section lift coefficient 
d e nozzle -exit diameter 

Kj slope of line-source or line-sink distribution of interval I 

l c length of chordwise vortex -line segment on wing 

fcore length of core of jet exhaust plume 

M Mach number 

N number of points or number of elemental vortex -lattice panels 

p static pressure 

Pt stagnation pressure 

Q strength of line-source or line-sink distribution 

q dynamic pressure 

S re f reference area (wing planform area as shown in fig. 1) 

s semiwidth of horseshoe vortex 

U longitudinal velocity 

u i> v i’ w i circulation-induced perturbation -velocity components in positive X-, Y-, and 
Z -directions 

u jet’ v j e t’Wj e t jet-induced perturbation-velocity components in positive X-, Y-, and 
Z -directions 

V e effective velocity ratio, 
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radial velocity in jet axis system induced by line sink- source distribution 


V r 


X,Y,Z orthogonal right-handed primary Cartesian coordinate system with origin at 
wing apex (see fig. 1). Positive X-direction is forward, positive Y-direc- 
tion is toward right wing tip, and positive Z-direction is downward. 
XY-plane parallel to wing chord plane 

Xj e t,Rj e t 3 et cylindrical coordinate system with origin at nozzle exit (fig. 3); positive 
X-direction is in direction of jet exhaust 


x,r cylindrical coordinates associated with jet (fig. 3), where a subscript indi- 

cates a specific value of each coordinate 


x c /4 longitudinal location in primary Cartesian coordinate system of midspan 

point of quarter -chord line of elemental panel (fig. 2) 


x 3c/4 


longitudinal location in primary Cartesian coordinate system of midspan 
point of three-quarter-chord line of elemental panel (fig. 2) 


x e ,y e ,z e primary Cartesian coordinates of nozzle exit (fig. 1) 
x^ e primary Cartesian X-coordinate of wing leading edge at y e (fig. 1) 

y C p primary Cartesian Y-coordinate of control point on elemental panel (fig. 2) 

a angle of attack, degrees 


r 

r 

C 

v 

e 

A 


circulation strength 
ratio of specific heats 
dummy variable , of integration 

fraction of wing semispan from wing-root chord in Y-direction 
camber angle of wing at control point, degrees 

sweep angle of quarter-chord line of elemental vortex -lattice panel (fig. 2), 
degrees 
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A c /4 

sweep angle of wing quarter-chord line, degrees 

Me 

wing leading- edge sweep angle 

X 

taper ratio, c^/c r 

P 

density 

Subscripts: 


c 

chordwise 

e 

exit 

r 

root 

s 

spanwise 

t 

tip 

total 

total 

oo 

free stream 


METHOD OF ANALYSIS 

The algorithm used in the present study to calculate the effects of exhaust jets 
located above a swept tapered wing is based on a vortex- lattice representation of the 
wing lifting surface and a line source- sink distribution to represent the entrainment and 
blockage effects of the exhaust jets. The method, which has been programed for a digital 
computer, assumes that the flow is steady, irrotational, and incompressible and that the 
flow external to the jet exhaust is inviscid. It is also assumed that the jet is not deflected, 
that the cross-sectional shape of the jet is not distorted by the wing flow field or by the 
free-stream angle of attack, and that the jets do not interact or wash the wings. No 
effects of the wing flow field on the entrainment of the jet were considered. Shown in 
figure 1 is a typical wing planform and engine arrangement for which the effects of 
blowing jets over the wing can be calculated by the present method. The present com- 
puter program does not allow discontinuous wing leading or trailing edges, but these 
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planforms are not excluded from the present method. Note also that the effects of the 
nacelles on the wing are neglected with the present method. 

The vortex- lattice method used to represent the wing lifting surface is based on 
references 13 and 14. For the present simplified planforms, each half of the wing is 
subdivided in both the spanwise and chordwise directions into 50 or less elemental areas. 
In the present computer program the number of spanwise and chordwise panels is arbi- 
trary; however, unless otherwise specified, all calculations for the present study were 
made with 5 equally spaced chordwise and 10 equally spaced spanwise divisions on 
each half of the wing. Each elemental area is represented by a horseshoe vortex with 
the bound portion lying along the local quarter-chord line of the element. (See fig. 2.) 
The trailing vortices lie streamwise along the inboard and outboard edges of each panel 
on the wing- chord plane. Aft of the wing trailing edge, the trailing vortices continue in 
the chordwise direction; that is, it is assumed that the trailing vortices are not turned bj 
the wing downwash or by the free-stream velocity. Note that this is not the conventional 
method of locating the trailing vortices parallel to the free-stream velocity, as used in 
reference 14, but is similar to the formulation of reference 13. The boundary condition 
that the flow be tangential to the wing surface is satisfied for each element at a control 
point on the lateral midpoint of the local three-quarter-chord line of the element. 

The induced effects of the jet exhaust were simulated with a distribution of line 
sinks and sources located on the longitudinal axis of the jet. The axis of the jet was 
divided into a number of segments over which the sink or source strength was assumed 
to be linear. This continuous line sink- source distribution is equivalent to a series of 
triangular elemental distributions as illustrated in figure 3. The radial velocity induced 
by the jet at a point in the surrounding flow field can therefore be obtained by summing 
the induced velocities resulting from each triangular elemental line source or sink. The 
radial velocity is therefore given by 

V (x _ Kir r x 2 ( x 2 - ^ ^ N y Kir f( x l- x l-l) r x bl ( x hl - g)dg 

1 ’ ‘ 41 \ ftx - 0* ♦ r*] 3/2 k 4 ’ ifl+i - X l) Jx I [(x - o 2 ♦ rf /2 

r x i (e- x i-i) d c %r «t N (; - * H )d; 

X I-1 [(x - f) 2 + r 2 ] 3 '' 2 / 4,1 X N-1 [(x - ri 2 + r^j 3,/2 

The slope Kj (I * 1,N) of the line source or sink of each elemental distribution was 
adjusted to give the inflow velocity on the boundary of the jet predicted by the method of 
reference 12. In the present analysis only those interference effects resulting from the 
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jet radial perturbation velocities are considered; any perturbation velocities induced 
parallel to the jet axis are neglected. In essence, then, the interference effects of the 
jet on the wing are attributed to an upwash distribution along the wing lifting surface. 

In order for the circulation strength of each elemental horseshoe vortex on the 
wing to be determined, the following equation for the circulation-induced downwash 
velocity must be satisfied at each control point on the wing: 

ZL = sm( tt H- . get _aet tan6 
cos 9 U„ U„ 

A sketch illustrating this boundary condition at each control point is shown in figure 4. 

It is apparent from this sketch how the induced effects of the jet modify the boundary- 
condition equation. Note that in the present formulation, the uj e t perturbation velocity 
arises only when the axis of the jet is not parallel to the wing- chord plane. 

After the circulation strengths for the horseshoe vortices are determined, the 
forces acting on the wing can be determined. From reference 13, the forces acting on 
each spanwise segment of vortex filament in coefficient form are 


and 


C A ,s = 



U 


oo 



+ 



- sin a) 



C 


2T, 

N > s 



+ 


tan A - 


Jh 


Vt 


+ cos a 



For each chordwise segment of vortex filament, the force coefficients are 
C A,c = 0 

and 


c = 2 rc f v i . V i e A 
N ’ c U TO iUoo ujs r ef 

where T c is the net circulation strength resulting from the individual circulations of 
each trailing horseshoe-vortex leg forming that segment. The total normal force 
Cn, total anc * axial force Ca, total ac ^ing on the wing are obtained by summing the 
normal and axial forces acting on all segments. The lift and induced-drag coefficients 
are then obtained from 
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g L = G N, total cos “ - C A, total sin « 

and 

C D,i = C A, total cos “ + C N, total sin « 

This solution for the induced drag is equivalent to the near-field solution of refer- 
ence 14. The integration of the wing span load distribution (the far-field method) to 
obtain the induced drag is not applicable to the present problem because the induced 
upwash field generated by the jet exhaust violates the assumption inherent in the far- 
field method. (See refs. 15 and 16.) 


DISCUSSION 

Comparison of Present Method With Other Theoretical Methods 

and With Experiment 

To verify the vortex- lattice algorithm used to represent the wing in the present 
method, the computed lift and drag characteristics of an aspect- ratio-8 straight wing 
without jets have been compared to the predictions of reference 14. (See fig. 5.) The 
present computer algorithm gives essentially the same results as reference 12 for each 
of the cases considered: a flat wing, a wing with Clark Y camber distribution, and a 
twisted wing with no camber. (The twisted wing had a linear distribution of twist which 
varied from -0.25° at the wing root to 4.75° at the wing tip.) Note that the near -field 
induced-drag results of reference 14 are compared with the present method. 

Experimental investigations of the effects of blowing jets over a wing are very 
scarce for the case in which the jet is sufficiently high to insure that it does not wash 
the wing. One source of such data is reference 17, where one jet is centered above a 
wing. The comparison of the present method with these experimental results, shown in 
figure 6, indicates very good agreement at vertical locations of the jet exhaust greater 
than 1.5 nozzle-exit diameters above the wing. In general, the effects of the ratio of jet 
velocity to free- stream velocity and of the vertical distance of the jet above the wing on 
the interference lift are well predicted. As the jet approaches the wing, however, the 
experimental results show a significant increase in interference lift over the predicted 
values. This discrepancy results, in part, from the interference lift associated with 
Coanda turning of the jet when the jet washes the upper surface of the wing. The effect 
of Coanda turning of the jet exhaust is, of course, neglected in the present method. 

A comparison of the present method with the experimental data of reference 9 (for 
two air jets blowing over a wing having a 50° swept leading edge with an aspect ratio of 3 
and a taper ratio of 0.3) is shown in figure 7. These experimental data were obtained at 
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ratios of nozzle -exit total pressure to free -stream static pressure from approximately 
1 to 6 and at free -stream Mach numbers of 0.4, 0.6, and 0.7. The theory of reference 10 
was derived for an incompressible jet exhausting into an incompressible external stream 
with no density difference between the two streams. The variable-density case can be 
related to the incompressible constant-density solution of reference 12 with an effective 
velocity ratio (ref. 18) as defined in the following equation: 



For an air jet with y g = = 1.4, 

M„(l + 0.2M e 2 f /4 

Ve = M7 

As shown in figures 7(a) and 7(b), the effective velocity ratio reasonably correlates 
the interference lift and the interference induced drag resulting from blowing two air 
jets over the wing having a 50° swept leading edge (ref. 9) at Mach numbers from 0.4 to 
0.7. The predictions of the present method are in good agreement with the experimental 
results. These results indicate that the present method can be extended by use of the 
effective-velocity-ratio concept to the case in which the jet and free stream are com- 
pressible and are at different densities. (Note that the exit-momentum coefficient of the 
jet is inversely proportional to the square of the effective velocity ratio and should pro- 
vide a similar correlation of the data.) 

A further application of the present method using the effective -velocity-ratio con- 
cept is shown in figure 7(c) for the 50° swept wing of reference 9. Here the following 
equations were used to obtain the predicted lift and induced drag: 




predicted 


L,jet offj measurec j 



predicted 


at a given a and 

(^T),i)predicted (^D,i)jet offj measure( j + i)predicted 

at a given Cl. The predicted lift and drag characteristics using the present method to 
calculate the effects of the jets with these equations are in good agreement with the exper- 
imental results. Of particular interest is the agreement of the predictions of the present 
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method with the experimentally indicated improvement in induced drag (not necessarily 
in total drag) due to blowing two jets over the wing. 

Effect of Vortex- Lattice Arrangement 

, Shown in figure 8 is the effect of the number of elemental chordwise panels on the 
predicted lift and induced drag of an aspect-ratio-8, taper-ratio-0.3 wing having the 
quarter chord swept 0° with and without two jets blowing over the wing. It is apparent 
that there is essentially no effect of varying the number of chordwise elemental panels 
on the predicted lift and induced drag of the wing or on the increments in lift and dr'ag 
due to jet blowing. Varying the number of spanwise elemental panels does result in a 
change in the predicted lift and induced drag for the wing with and without jet blowing. 

(See fig. 9.) It appears, however, that the interference effect due to jet blowing is not a 
strong function of the number of elemental vortex- lattice panels. Therefore, within the 
limits of the present computer program, the predicted effects of jet blowing are essen- 
tially unaffected by vortex-lattice arrangement. For that reason, 5 chordwise and 
10 spanwise elemental panels have been used for all calculations of the present paper. 

Effect of Velocity Ratio 

The predicted effects of effective velocity ratio on the interference lift and inter- 
ference induced drag due to jet blowing are shown in figures 6, 7, and 10 for an aspect- 
ratio-2 straight untapered wing, for an aspect- ratio- 3, taper-ratio-0.3 wing with a 50° 
swept leading edge, and for an aspect- ratio-8, taper- ratio-0. 3 wing with a 30° swept 
quarter chord, respectively. The present method predicts an increase in the interference 
lift coefficient with increasing jet blowing (that is, increasing 
interference effect is essentially unaffected by an angle of attack of up to 10°. (See 
fig. 10(c).) The present method also predicts a decrease in induced drag with an increas- 
ing ratio of jet velocity to free-stream velocity l/V e . This favorable increment in 
interference-induced drag increases with the lift coefficient. (See fig. 10(c).) In 
figure 10(b) are shown typical calculations of the effects of effective velocity ratio on the 
wing span load distribution and on the section lift coefficient. Increased jet blowing 
results in an increased distortion of the span load distribution. Note that this increased 
distortion does not imply increased induced'drag since the integration of the span load 
distribution to obtain induced drag is only valid for planar wings without external inter- 
ference effects. As would be expected, the span load does increase in the vicinity of the 
jet. Increasing jet blowing also increases the section lift coefficient over the entire 
wing span. 


l/V e ). This favorable 
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Effect of Jet Location 


The predicted effects of jet location on the interference lift and interference induced 
drag caused by one or two jets blowing over various wing configurations are shown in 
figures 6 and 11 to 14. Increasing the vertical distance of the jets above the wing (figs. 6 
and 11) causes a decrease in the interference lift and a reduction in the favorable inter- 
ference induced drag. Moving the nozzle exit ahead of the wing leading edge (fig. 12) 
causes a small increase in the interference lift and a small increase in the favorable 
induced-drag increment. Moving the nozzle exit aft of the wing leading edge causes a 
decrease in the interference lift and a substantial reduction in the favorable induced-drag 
increment. The calculations shown in figures 13 and 14 indicate that the closer the jets 
are located to the wing center line the more favorable the interference effects caused by 
blowing two jets over a wing. This is indeed a fortuitous result, since reference 4 indi- 
cates that for a two- engine configuration the engine should be located inboard to mini- 
mize the engine-out effect on the rolling moment associated with upper surface blowing 
configurations during take-off and landing. 

CONCLUDING REMARKS 

A procedure has been developed for calculating the effects of blowing two jets over 
a swept tapered wing at low subsonic speeds. The algorithm used is based on a vortex- 
lattice representation of the wing lifting surface and a line sink- source distribution to 
simulate the effects of the jet exhaust on the wing lift and drag. The method is limited 
to those cases in which the jet exhaust does not intersect or wash the wing. The pre- 
dictions of this relatively simple procedure are in remarkably good agreement with 
experimentally measured interference lift and interference induced drag. 

The results of the analytical study have indicated that substantial increases in 
wing lift and reductions in wing induced drag can occur when jets are blown over a wing. 
The magnitude of these effects increases with an increasing ratio of jet-exit velocity to 
free-stream velocity. The results also indicate that the interference effects are the 
most favorable when the nozzle exits are located ahead of or at the wing leading edge 
and inboard near the wing center line. The present method does not correctly predict 
the effects of the vertical location of the jets at jet center-line heights of less than 
approximately 1.5 jet-exit diameters; however, as the vertical location of the jets is 
reduced to this height, the favorable effects of blowing jets over a wing increase. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 21, 1974. 
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gure 3.- Illustration of source-sink distribution used to represent effects of jet exhaust. 
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line) 

Boundary condition: 

UcpSin (a + 0) = (wj + w- et ) cos 0 + u. gt sin 0 

Figure 4.- Sketch illustrating boundary conditions at each control point 

on a wing with camber. 
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(a) Flat wing. 


Figure 5.- Comparison of present method with that of reference 14 for an 
aspect-ratio-8, A. = 1.0 straight wing (jet off). (Near -field solution of 
ref. 14 for induced drag is shown.) 
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(b) Wing with Clark Y camber distribution. 
Figure 5.- Continued. 
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(c) Wing with twist distribution. 
Figure 5.- Concluded. 





Figure 6.- Comparision of present method with experimental data of reference 17 


for blowing one jet over an aspect-ratio-2, A. = 1 straight wing. A e/ /S re j = 0.192, 
( x e - x Ze)A> = 1*0, V e = °-°> c e/ d e = 6 - 250 > and a = °°- 
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(a) Interference lift coefficient. 

Figure 7.- Comparison of the present method with experimental data of reference 9 
for two jets blowing over a wing with aspect ratio 3, X = 0.3, and Aj e = 50°. 
Ae/Sref = ^*0078, -/ x e - x Ze)/^e = 2.59, rj e = 0.46, and z e /d e = -1.5. 














MqO 


O 0.4 "j 

□ .6 f Experiment 



2 3 4 5 


'/ v e 

(b) Interference induced drag coefficient. 
Figure 7.- Continued. 
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Figure 8.- Effect of number of chordwise elemental vortex -lattice panels when 
N s = 10 on the calculated lift and induced drag of a wing with aspect ratio 8, 
\ = 0.3, and A ^4 = 0° and with two jets having A e /S re f = 0.025 located at 
fx e - xj e Vcl e = 0.0, r] e = 0.3, z e /d e = -1.5, and c e /d e = 2.408. 
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Figure 9.- Effect of the number of spanwise elemental vortex -lattice panels when N c = 5 on the calculated lift and 
induced drag of a wing with aspect ratio 8, X = 0.3, and A c / 4 = 0° and with two jets having A e /S ref = 0.025 
located at (x e - xj e )/d e = 0.0, rj e = 0.3, z e /de = -1.5, and e e /d e = 2.408. 
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(a) Effect on C L and Cj)^. 

Figure 10.- Calculated effect of effective velocity ratio of two jets on aerodynamic char- 
acteristics of a wing with aspect ratio 8, \ = 0.3, and A c /4 = 30°, A e /S re f = 0.0i25, 
(x e - He) Me = 0.0, r] e = 0.25, z 6/ /d e = -1.401, and c G/ /d e = 3.557. 
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(b) Effect on span load distribution and section lift coefficient 

at a = 2°. 


Figure 10.- Continued. 
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(a) Effect on and Cj^. 

Figure 11.- Calculated effect of vertical location of two jets on aerodynamic character- 
istics of an aspect-ratio-8, A. = 1.0 straight wing. A e/ /S re f = 0.0125, V e = 0.20, 
rj e = 0.25, ^x e - l e = 0.0, and c e /d e = 2.803. 
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(b) Effect on span load distribution and section lift coefficient 

at a = 2°. 


Figure 11.- Continued. 
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(c) Effect on ACl and ACj-^. 
Figure 11.- Concluded. 
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C L 

(a) Effect on Cj^ and Cj^j. 

Figure 12.- Calculated effect of longitudinal location of two jets on aerodynamic 
characteristics of an aspect-ratio-8, \ = 1.0 straight wing. A e /S re f = 0.0125, 
V e = 0.20, t] e = 0.25, z e /d e = -1.401, and c e /de = 2.803. 
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(b) Effect on span load distribution and section lift coefficient 

at a = 2°. 


Figure 12.- Continued. 
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(a) Effect on and Cp 

Figure 13.- Calculated effect of lateral locat*on of two jets on aerodynamic char- 
acteristics of an aspect-ratio-8, \ =* 1.0 straight wing. A e /S re f = 0.0125, 

V e = 0.20, ^x e - dg = 0.0, z e /d e = -1.401, and c e /d e = 2 -803. 
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(b) Effect on span load distribution and section lift coefficient 

at a = 2°. 


Figure 13.- Continued. 
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(c) Effect on ACl and ACp^. 
Figure 13.- Concluded. 
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(a) Effect on Cl and Cp j. 

Figure 14.- Calculated effect of lateral location of two jets on aerodynamic 
characteristics of a wing with aspect ratio 8, X = 0.3, and A c /4 = 30°. 
A e /S re f = 0.0125, V e = 0.20, ^x e - xj e ^/d e = 0.0, and z e /d e = -1.401. 
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(b) Effect on span load distribution and section lift coefficient 

at a = 2°. 

Figure 14.- Continued. 
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(c) Effect on ACl and AC 
Figure 14.- Concluded. 
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